Elucidation of the anomalous A = 9 isospin quartet behaviour 
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Recent high-precision mass measurements of 9 Li and 9 Be, performed with the TITAN Penning 
trap at the TRIUMF ISAC facility, are analyzed in light of state-of-the-art shell model calculations. 
We find an explanation for the anomalous Isobaric Mass Multiplet Equation (IMME) behaviour for 
the two A = 9 quartets. The presence of a cubic d — 6.3(17) keV term for the J 71 = 3/2~ quartet and 
the vanishing cubic term for the excited J 71 = 1/2" multiplet depend upon the presence of a nearby 
T — 1/2 state in 9 B and 9 Be that induces isospin mixing. This is contrary to previous hypotheses 
involving purely Coulomb and charge-dependent effects. T — 1/2 states have been observed near 
the calculated energy, above the T = 3/2 state. However an experimental confirmation of their J n 
is needed. 
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Atomic nuclei are described by their binding energy 
and three quantum numbers: the total angular momen- 
tum J, parity 7r, and isospin T. This framework allows 
one to identify, each of the ~ 3000 observed nuclei [l[ un- 
ambiguously. The isospin quantity is analogous to spin 
and was first introduced by Heisenberg to describe 
the charge-independence of the nuclear force. Within the 
isospin formalism, neutrons (n) and protons (p) are nu- 
cleons of isospin T — 1/2 but distinguished by different 
z-projections T 2 (n) = 1/2 and T z (p) = -1/2 @,[3|. Nuclei 
with the same mass number A, total angular momentum 
and parity form multiplets where the individual members 
have a projection T z = (N — Z)/2. Assuming isospin is 
a good quantum number, members of an isobaric mul- 
tiplet have identical properties. However, Weinberg and 
Treiman [1] noted that the mass excess A (which is a 
measure of the nuclear binding energy and defined as the 
difference between the atomic mass and the atomic mass 
number) of such nuclides were not identical, but were 
rather laying along a parabola: 

A(A, T, T z ) = a{A, T) + b{A, T)T Z + c{A, T)T 2 Z (1) 

where a, 6, c are coefficients that depend on all quan- 
tum numbers except T z . This so-called isobaric multiplet 
mass equation (IMME) has proven to be a powerful tool 
to predict unknown masses. For instance, it is used to 
obtain masses of nuclei along the rapid proton capture 
path, where most of the masses are not well known |5j or 
to provide detailed mass values, which are experimentally 
inaccessible due to half-life and productions constraints 
Q. Recently, the precise mass measurement of 12 Be Q 
using the TITAN (TRIUMF Ion Traps for Atomic and 
Nuclear science) Penning trap mass spectrometer 0, Q 
has been used as a solid anchor point together with the 
IMME to address the ambiguous spin assignment of T = 
2 states in 12 C and 12 Be. 



Several tests of the IMME were performed and for most 
cases, it has followed the original quadratic behaviour 
[irj| . However, in some cases, large deviations were dis- 
covered and the incorporation of cubic d(A, T)T Z and/or 
quartic e(A,T)T z [HI, EH terms was considered. The 
largest documented breakdowns include the A = 9 [IH 
and A = 33 quartets and the A = 8 [l6| and A = 
32 [l7, 18 1 quintets. The unveiling of the non-quadratic 
behaviour of the A = 32 and 33 multiplets was only pos- 
sible due to the precise and accurate mass measurement 
of some of its members, at the 5m/ m ~ 10~ 7 level, us- 
ing Penning traps [l9j]. Because of their lighter mass, 
the A = 8 and 9 multiplet breakdowns where discovered 
sooner, from less precise reaction Q- value mass detcrmi- 
*2ol . 21 1 . More recently, mass measurements of 
22j and 8 He [3] showed the need for a larger cubic 
d = 11.1(2.3) keV term in the A = 8 quintet [3]. In 
addition, a new IMME evaluation of the A — 9 quartet 
became possible using recent lithium |2J] and beryllium 
25j mass measurements. The A = 9 isobars are of par- 
ticular interest because it is the first and ligh test chain 
presenting two different IMME quartets The 



nations 
8 C 



ground state quartet [21| strongly departs from quadra- 
ture with a quadratic fit x 2 of 10.2 and a cubic coefficient 
d = 5.5(18) keV [HI. Several mechanisms to explain this 
departure have been proposed including the Coulomb- 
dependent Thomas-Ehrman shift |27| arising from the 
small binding energy of the last proton in 9 C and non- 
Coulomb charge-dependant forces 12|. However, to date, 
the total contribution from these mechanisms was insuf- 
ficient to explain the observed cubic term 
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thermore, as the strength of these effects would increase 
with a decreasing proton separation energy in 9 C, the 
excited state quartet should depart more strongly from 
quadrature [12j. However, the excited state quartet in A 
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FIG. 1: (Color on-line) New mass excess A(new) values pre- 
sented in tablcUcompared to the 1995 atomic mass evaluation 
(AME1995) H mass excess values A(AME1995) used in the 
most recent IMME review [ljj. The new 9 Li and 9 Be A were 
measured at TITAN, while the 9 B and 9 C A involve respec- 
tively a 9 Be(p,n) 9 B and 7 Be( 3 He,n) 9 C reaction Q-value. 



= 9 shows a good agreement with quadrature (x 2 = 1-0) 
and a nearly vanishing cubic term d = 3.5(34) keV [Ioj [. 
The mixing of the T = 3/2 and T = 1/2 states for the T z 
— ±1/2 members has been proposed [ll[, but the large 
width of the known resonances met with skepticism [28[ 
for the strength of this mechanism. This Letter explains 
the long-standing enigmatic behaviour of the two A = 9 
quartets. 



TABLE I: New ground state mass excesses Acs. and excita- 
tion energies [l[ of the T — 3/2 states used for the calculation 
of the r = 3/2" and 1/2" A = 9 quartets. The 9 B ,F = 
1/2" excitation energy labeled with (*) is from [Til ]. 



T 2 



Acs. 
(keV) 



(keV) 



(keV) 



J Li 3/2 24 954.91(20) 2 691(5) 

'Be 1/2 11 348.391(93) 14 392.2(18) 16 977.1(5) 
5 B -1/2 12 416.4(10) 14 655.0(25) 17 076(4) 

16 990(30) * 

J C -3/2 28 909.5(21) 2 218(11) 



The most recent IMME review (1998), used the ground 
state masses from the 1995 atomic mass evaluation 
(AME1995) 0. In this evaluation, all ground state 
masses for the A = 9 quartets where based on reaction Q- 
value measurements from transfer reactions. Since then, 
the masses of 9 Li 24| and 9 Be 
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rectly using the TITAN Penning trap mass spectrom- 
eter of the ISAC facility at TRIUMF. Penning traps 
have been established as the most precise and reliable 



devices for mass measurements [19(, which is a key pre- 
requisite for this IMME study. The TITAN system has 
been established with mass measurements of halo nuclei 
as well as with measurements of short-lived, 
highly-charged ions j3l[. The precision and accuracy 
of the TITAN Penning trap has been studied in detail 



The new masses are presented in table HI The large 
improvement in precision and change in the mass value 
compared to AME1995 is displayed in Fig. [TJ The TI- 
TAN 9 Be mass excess also affects 9 B, which is derived 
from the Q-value of a 9 Be(p,n) 9 B reaction [33, |35|. The 
4.2 keV decrease in the 9 C mass excess, currently de- 
rived from a 7 Be( 3 He,n) 9 C Q- value, originates from the 
change in the 7 Be mass evaluation [34j,|35j . The 7 Be mass 
value in the 2003 Atomic Mass Evaluation (AME03) was 
based solely on the Q-value of the reaction 7 Li(p,n) 7 Be 
[35| which was not the case for the 1995 evaluation f34j ]. 
The 9 C mass presented in table U is also based on this 
reaction, but uses the precise 7 Li mass value measured 
using the SMILETRAP Penning trap spectrometer (3(|. 
All the excitation energies presented in table U come from 
the Evaluated Nuclear Structure Data File (ENSDF) Q, 
except for the 9 B J T = 1/2" state. For this state we 
include a recent 3<r lower value of 16 990(30) keV [22| 
together with the 17 076(4) keV value from ENSDF. In 
the new measurement, the angular correlation between 
the first proton emission of 9 B and subsequent double a 
decay strongly suggest the J 71 ' = 1/2" assignment for this 
state. 



TABLE II: IMME coefficients (equation (fl|) using the mass 
excesses calculated from the values in table H] Also given are 
the x 2 an d P-value of the various fits. The (*) value used the 
9 B J 71 " = 1/2" excitation energy from [22l |. 
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3/2" 


26337.5(17) 


-1318.8(7) 264.7(9) 14.6 


1X10" 4 


1/2" 


28847.1(15) 


-1163.7(29) 241.2(24) 1.2 


3X10" 1 


1/2" 


* 28845.0(19) 


-1159.1(41) 240.1(25) 6.9 
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The updated IMME coefficients for the fits, including 
quadratic and cubic terms are shown in Table ITT1 together 
with their x 2 ■ As the various fits involve only one degree 
of freedom, we also give the corresponding P-value for the 
X 2 ■ The new mass excesses result in a 40% increase in the 
2 for the J 71 ' = 3/2" quartet quadratic fit compared to 



X 

the previous value [lfj. This very large x 2 is associated 
with a small P-value, which reflects the probability to ob- 
tain such a x 2 by accident. Performing a cubic fit results 
in an enhanced d coefficient of 6.3(17) keV. The excited 
J 71 = 1/2" quartet, on the other hand, shows a good 
agreement with quadrature leading to ay 2 = 1.2 when 
using the ENSDF 9 B excitation energy [1(. However, if 
the 16 990(30) keV value for the 9 B r = 1/2" excita- 
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TABLE III: Experimental and theoretical d, 61 and 63 coeffi- 
cients for the two A = 9. 
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FIG. 2: (Color on-line) Experimental and theoretical energies 
for negative parity levels of 9 B. The total angular momenta 
J are indicated by the length of the line. The experimental 
levels for which J is uncertain are indicated by the black dots. 
The states indicated by a thick red line ending with a circle 
are the T = 3/2 states of interest. The theoretical state in- 
dicated by a thick blue line ending with a cross is the T = 
1/2 state that isospin mixes with the lower T — 3/2 state. 
The experimental state in blue is a possible candidate for the 
associated T = 1/2, J n = 3/2~ level. 



able data set for fitting. The non-zero cubic term for the 
J n = 3/2" , T = 3/2 level is dominated by isospin mixing 
with a higher J w = 3/2", T = 1/2 state. Calculations 
using the CKI Hamiltonian predicts this level to be about 
0.5 McV higher, while the PJT Hamiltonian predicts this 
level to be about 0.8 MeV higher (right-hand side Fig. [2]). 
Experimentally, aT = 1/2 level with unknown total an- 
gular momentum, parity and width is observed very close 
to the T = 3/2 level expected from PJT, at 0.6 MeV 
(left-hand side of Fig. [2]) [l|. This isospin mixing can 
be understood as a second-order perturbation in which 
these two levels repel each other. As a result, the T = 
3/2 level is pushed down by an amount proportional to 
the expectation value of the isospin non-conserving in- 
teraction divided by the energy difference AE between 
the two states: |< Vinc >| 2 / AE. The difference in the 
predicted d coefficient from the CKI and PJT Hamilto- 
nian is mainly due to the different energy denominator. 
If the mixing originates from the observed T = 1/2 state 
suggested above, then the PJT results would be in good 
agreement with experiment for both AE and d. In con- 
trast, calculations using either Hamiltonian do not pre- 
dict a J n = 1/2- , T =1/2 state near the quartet J* = 
1/2", T = 3/2 state, leading to a nearly zero d coeffi- 
cient. This is confirmed experimentally as no nearby T 





6l(3/2-) 63(3/2-) rJ(3/2") 
(keV) (keV) (keV) 


6i(l/2-) b 3 (l/2-) d(l/2~) 
(keV) (keV) (keV) 


exp 
PJT 
CKI 


-1330.8(32) -1318.2(7) 6.3(17) 
-1334 -1321 6(2) 
-1364 -1342 11(2) 


-1161(4) -1167(4) 3.2(2.9) 
-1258 -1260 -1(2) 
-1290 -1290 0(2) 



tion energy from 22| is used, the excited state quartet 
also departs from quadrature with & x 2 — 6.9 and a larger 
cubic term of d = -40(15) keV. A more precise confirma- 
tion of this recent measurement is desired, to resolve this 
discrepancy for the excited quartet. 

In order to see how the large cubic term of the ground 
state quartet arises, we define two new b coefficients, b\ 
and 63: 

h = A(A, T z = 1/2) - A(A, T z = -1/2) (2) 
b 3 = [A(A, T z = 3/2) - A(A, T z = -3/2)]/3. (3) 

The b and d coefficients can be written in terms of 61 
and & 3 as d = [61 - 63] /2 and b = [%i - 63] /8. From 
these expressions, it can be seen that for 63 = b±, the 
other coefficients become d = and b = 63 = b\. With 
this change, the physics involving only the T z = ±1/2 
states is now explicitly contained in the single b\ term. 
The experimental and calculated results for these coef- 
ficients are presented in table IIIII The results were cal- 
culated using a shell model Hamiltonian comprising two 
parts: isospin conserving and isospin non-conserving [37J 
component. The first part uses either the Cohen-Kurath 
two-body matrix clement (6-16)2BME (CKI) [H] or the 
updated version of this, the so-called p-shell potential 
model fit (PJT) [3<|. The PJT Hamiltonian is given in 
[io| . Both of these are based on fits of the two isospin- 
conscrving single-particle energies and fifteen two-body 
matrix elements adjusted to binding energies and excita- 
tion energies for the p-shell nuclei from A — 6 to 16. The 
PJT Hamiltonian takes into account more data than was 
available for the CKI. The isospin non-conserving Hamil- 
tonian includes a Coulomb, charge-asymmetric isovector 
and charge-dependent isotensor strong interaction. The 
free parameters of these interactions were obtained by 
fitting the Hamiltonian to the p-shell nuclei experimen- 
tal data for the b and c coefficients of the IMME. The 
calculated energy levels, total angular momentum and 
parity using the PJT Hamiltonian are compared to the 
experimental values in Fig. [2] The following discussion 
will focus on the J x = 1/2" and 3/2", T = 3/2 states 
indicated by thick red lines in Fig. [2] 

Tabic ITTTl shows that the theoretical results for the d co- 
efficient, calculated to one keV precision, are in excellent 
agreement with experiment. The better agreement for 
the PJT potential could be attributed to its larger avail- 
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= 1/2 states have been observed 

The interference between the Coulomb, charge- 
asymmetric isovector and charge-dependent isotensor 
part of the isospin non-conserving interaction is the crit- 
ical ingredient that gives rise to the level repulsion stem- 
ming from the isospin mixing. If the Coulomb term is put 
to zero and either the isotensor or isovector part of the in- 
teractions is used, then the corresponding matrix element 
is the same for both 9 Be and 9 B, resulting in an equal 
downward push of their respective T = 3/2 states and a d 
= term. In detail, the calculation decomposes these ma- 
trix elements into a sum of Coulomb, charge-asymmetric 
isovector and charge-dependent isotensor contributions 
resulting in |< V mc >| = I 54 + 44 + 18 | = 116 keV 
and |< V mc >| = | - 19 - 44 + 18 | = 45 kcV for 9 B 
and 9 Be respectively. The square of these values divided 
by the energy difference between the T = 3/2 and nearby 
T = 1/2 state yield downward shifts of about 14 keV in 
9 B and 3 keV in 9 Be. Then, the d coefficient is equal to 
half the difference between the two downward shifts. 

The isospin non-conserving matrix elements theoreti- 
cal uncertainties are related to the model and uncertain- 
ties used for fitting the b and c coefficients as given in Ta- 
ble 2 of [33] • The Coulomb term is well determined and 
is within 5% of its expected strength. The c-coefheient 
is sensitive to the charge-dependent isotensor term fl3j ]. 
Its value and error are given in Table 2 of [37[ as Sq 



-0.017(5), leading to an error of about 1 keV in the d- 
coefficient. The charge-asymmetric term was determined 
from a fit to the fe-coefficients and given in Table 2 of [37j 
as Sq 1 ^ = -0.042(11), leading to an error of about 1 keV 
in the d-coefficient. Both of these contributions yield a 
total theoretical uncertainty of 2 kcV. In [37J these em- 
pirical coefficients were compared to simple model esti- 
mates based on the experimental nucleon-nucleon scat- 
tering lengths; S$ = -0.042(12) and = -0.006(4). Us- 
ing these estimates in place of the empirical values re- 
sults in d = 7 keV; essentially the same result although 
the details in terms of the division between the isovector 
and isotensor terms is rather different. In the future it 
would be desired that the isospin non-conserving contri- 
butions used for the | < Vinc > I should be reevaluated by 
rcnormalizing nucleon-nucleon interactions to the p-shcll 
model space with inclusion of three-body interactions. 

Previous work on the theoretical calculations of the 
IMME including the d coefficients is reviewed in [l3| . 
Calculations of the A = 9 J + = 3/2~ quartet was con- 
sidered in [l2[ where the Coulomb interaction was used 
in a schematic model, yielding d = 0.9 keV. Charge- 
dependent interactions resulted in an extra 1 keV. Most 
importantly, the results obtained in 



121 do not include 



In summary, combining high-precision 9 Li and 9 Be 
mass measurements performed using the TITAN Penning 
trap spectrometer, and shell model calculations using the 
most complete potential to date, we could fully explain 
for the first time the binding energy behaviour of the A = 
9 quartet. Contrary to previous hypothesis, higher-order 
charge dependant effects and Coulomb effect-induced ex- 
pansion of the wave function are not the leading mecha- 
nism driving a non-zero term in this quartet. It is created 
by the isospin mixing of the T = 3/2 level with an above 
T = 1/2 level. Using this approach, the long-standing 
anomalous difference in d-cocfficicnts for the two A = 9 
quartet is explained by the absence of a nearby T = 1/2 
level for the J 71 " = l/2~ quartet. Nevertheless, an ex- 
perimental confirmation of our predicted J T for the 15 
100(50) keV and 15 290(40) keV states in 9 Bc and 9 B as 
well as a measurement of their width would be of inter- 
est. Also, light should be shed on the excitation energy 
discrepancy for the 9 B J = l/2~ level. 

Finally, the presented interpretation for the creation of 
a non-zero d coefficient in the ground-state A = 9 quartet 
can be generalized to other isobar multiplets if the same 
conditions of a small energy separation between mixing 
levels and a non-negligible isospin non-conserving matrix 
element is met. This generalized theoretical description 
can also be used to explain other physical effects such as 
isospin-forbidden proton decays 4l| and the necessity of 
a correction term S c to properly determine the Cabibbo- 
Kobayashi-Maskawa matrix V u d mixing angle 18l. |42|. 
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the contribution from mixing with specific nearby states 
and they are thus not inconsistent with our result as our 
larger cubic term is primarily due to the mixing with a 
nearby state. 
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